Multi-Bloch Vector Representation of the Qutrit 
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An ability to describe quantum states directly by average values of measurement outcomes is 
provided by the Bloch vector. For an informationally complete set of measurements one can con- 
struct unique Bloch vector for any quantum state. However, not every Bloch vector corresponds to 
a quantum state. It seems that only for two dimensional quantum systems it is easy to distinguish 
proper Bloch vectors from improper ones, i.e. the ones corresponding to quantum states from the 
other ones. I propose an alternative approach to the problem in which more than one vector is 
used. In particular, I show that a state of the qutrit can be described by the three qubit-like Bloch 
vectors. 



I. INTRODUCTION 

A state of a discrete quantum system is described by 
a complex vector. If only partial information about the 
system is given, or if it is entangled with any other sys- 
tem, its state is represented by a density matrix. This 
well known description is always nonintuitive to anyone 
hearing for the first time of quantum mechanics because 
it takes place in the abstract Hilbert space instead of the 
spacctimc, in which all physical phenomena occurs. How- 
ever, some alternative representations of quantum states 
referring to classical intuition have been introduced. 

An important alternative description of a discrete 
quantum state is given by the Bloch vector which 
explicitly refers to experimental notions like expectation 
values of measurements which in principle can be per- 
formed in a laboratory. The state is represented by a real 
vector whose coordinates are simply expectation values 
of different measurements and whose dimension is given 
by a minimal number of measurements needed to obtain 
complete information about the system. Although widely 
used, the Bloch vector description is still under develop- 
ment 0, H, 0, 1 S 13 > because, although the idea seems 
simple, it is hard to provide one important feature. Every 
quantum state can be written as the Bloch vector, but 
not every vector corresponds to a quantum state. For 
some vectors the corresponding density matrix, though 
still of trace one, can have negative eigenvalues. 

In Hilbert space of dimension two the set of all Bloch 
vectors is confined in a real three dimensional ball. Such a 
nice symmetry has not been found for higher dimensions 
@, H H H !, 01 • The symmetry allows for an easy recog- 
nition of the proper Bloch vectors, i.e., the vectors corre- 
sponding to quantum states. However, there is also one 
more reason to look for well shaped symmetric bound- 
aries. The beauty of the Bloch ball is manifested by the 
fact that its symmetry not only represents nicely states of 
the qubit, but also gives very intuitive picture of the dy- 
namics of the system — every unitary operation is simply 
a rotation of the ball. Therefore, an additional goal of a 



generalization of the Bloch vector to higher dimensions 
is to give an intuitive representation of the dynamics of 
the underlying system. 

The Bloch vector of the qubit is well know to anyone 
working in the field of quantum information, optics or 
magnetism. On the other hand, it is obvious that ev- 
ery quantum system contains a two dimensional (qubit) 
subspace. In particular, the number of a distinct qubit 
subspaces of a d-level system is given by a binomial coef- 



ficient 



Each of these subspaces can be described by 



its own qubit Bloch vector. However, these Bloch vectors 
cannot be independent. 

An important extension of the qubit is the qutrit, 
the three level system. Despite apparent simplicity, the 
qutrit is complex enough to reveal bizarre properties 
of quantum theory, like contextuality ||. Moreover, it 
was shown that using qutrits, instead of qubits, one can 
perform much better in quantum information processing 
tasks, like for example the quantum key distribution [?|. 

The aim of this work is to find a set of constraints on 
Bloch vectors of qubit subspaces of the underlying qutrit. 
As a result, one obtains a new Bloch-like representation 
of the qutrit. Moreover, it is going to be shown that 
this representation is natural for spin 1 whose dynamics 
is intuitively fitted for this picture. The extension to 
higher dimensions is also discussed. 



II. OPERATOR BASIS 

The usual construction of the Bloch vector for an 
arbitrary qudit starts with a choice of an appropriate 
measurements which mathematically are represented by 
d 2 elements of an orthogonal hcrmitian operator ba- 
sis. The orthogonality is according to the trace product 
(A, B) = Tr{A'B}. Due to the normalization constraint, 
there are d 2 — l parameters characterizing the state of the 
qudit. Usually 0,1 11,10, 

all but one basis operator, 
which is just the identity divided by d, are chosen to be 
traceless so any density matrix can be written as 
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where Mi = Tr{Af } is the corresponding trace norm and 
at = Tr{pAi}. 

Let me introduce an operator basis for the qutrit, 
which is natural for a decomposition into distinct qubit 
subspaces. It consists of six traceless Gell-Mann matrices 
and three diagonal rank 1 matrices of trace 1. The 
first triple of Gell-Mann matrices is symmetric: 

A 12 = |1)(2| + |2)<1|, 
A 13 = |1)(3| + |3)(1|, 
A 23 = |2)(3| + |3)(2|. 

The second one is antisymmetric 

B 12 = -i|l)(2|+i|2)(l|, 
Bi 3 = -i|l)(3|+i|3)(l|, 
B 23 = -»|2)(3|+i|3)(2|, 

Both, symmetric and antisymmetric matrices have trace 
norm equal 2. The last triple is given by 

Oi = |l)(l|, 2 = |2)(2|, 3 = ]3)<3j. 

All nine operators form an orthogonal operator basis in 
which one can expand any 3x3 matrix. In particular, 
the density matrix can be written as 





(2) 



Throughout the paper it is assumed that for all subscripts 
* 7^ 3 7^ an d if for some subscript ij i > j, then it 
should be read as ji. The coefficients u>i, ctij and /3y are 
expectation values of the corresponding operators with 
respect to p. It is convenient to write Eq. @ in an 
explicit form 



ai2 +" 312 UJ 2 
Qi3+i/3i3 Q23-M0: 




The three 2x2 principal sub-matrices of the above ma- 
trix represent qubit subspaces. It is important to notice 
that for a positive semidefinite matrix all of its princi- 
pal sub-matrices are also positive semidefinite (see [loj). 
Thus, principal 2x2 sub-matrices of the qutrit can in- 
deed describe qubits, though unnormalized. The idea 
is to represent the qutrit by the three Bloch vectors of 
these qubits. As already mentioned, they cannot be inde- 
pendent, since qutrits are described by only 8 parameters 
and each qubit subspace is described by 3 + 1 parameters. 
An additional parameter is due to the lack of the nor- 
malization constraint, and the radius of the correspond- 
ing Bloch sphere can be less than 1. For the consistent 
picture one has to derive the full set of constraints on 
parameters u>i, cnj and fy. 



III. CONSTRAINTS ON BASIS EXPANSION 
COEFFICIENTS 

First of all, the trace of p has to equal 1, which gives 
the first constraint 



5> = i. 



(3) 
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Since coefficients tUi are the smallest principal sub- 
matrices of p, they have to be nonnegative. Moreover, 
due to relation each of them has to obey 



< u>k < 1. 



(4) 



There is nothing new about the above relation, because 
diagonal elements of density matrix correspond to prob- 
abilities of measurement outcomes in the computational 
basis. 

Next, let me find what conditions would provide the 
positive semi-definiteness of p. This can be done by 
studying principal minors with the help of the interlacing 
inequality [lCj . This inequality states that if A is an n x n 
hermitian matrix and A' is its (n — 1) x (n — 1) princi- 
pal sub-matrix, then the eigenvalues of A' fall exactly in 
between the eigenvalues of A 



Ai < a; < a 2 < ■ • • < a;_! < A, 



(5) 



The eigenvalues are indexed according to a nondecreas- 
ing order, i.e., if i < j then < Xj (the same holds 
for X'i). There exists a very nice criterion which states 
that A is positive definite if and only if all of its principal 
minors are positive (see [T(J). In fact, it is enough to 
consider only leading (upper left) principal minors. Un- 
fortunately, a similar criterion does not work for positive 
semidefinite matrices. However, below I show that for 
3x3 matrices of trace 1 the matrix p is indeed positive 
semidefinite if and only if all of its principal minors are 
nonnegative. 

If p is positive semidefinite, then by the interlacing 
inequality ([5]) all of its principal minors are nonnegative. 
To prove the converse, let me start with lxl principal 
sub-matrices, which are simply the coefficients u>{. Their 
non-negativity is already provided by (f4]). If some of 
them were negative, then by (J5j) p would have at least 
one negative eigenvalue. 

If a 2 x 2 principal minor is negative, then the corre- 
sponding 2x2 principal sub-matrix has one positive and 
one negative eigenvalue. This implies the negativity of 
at least one eigenvalue of p, thus no 2 x 2 principal mi- 
nor can be negative. If the same minor is positive, the 
interlacing inequality states that both eigenvalues of the 
corresponding sub-matrix have to be positive. Finally, if 
it is zero, one has to consider two cases. If both principal 
sub-matrices u>i and u>j are positive, then one eigenvalue 
is positive and one equals zero. On the other hand, if at 
least one u>i equals zero, then af • has to be zero too, 
but in this case a,j = fyj = 0, thus u>i and u>j are simply 
eigenvalues, which by definition are nonnegative. 
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The fact that 2x2 principal minors have to be non- 
negative can be formulated as three new constraints 



«7 



■0; 



(6) 



which resemble the well known constraint on the coeffi- 
cients of the qubit. Interestingly, the above constraints 
indicates the spherical symmetry on parameters aij and 

Since we already know that all principal sub-matrices 
of p, whose dimension is smaller than three, have to be 
positive semidefinite, what is left to study the positive 
semi-definiteness of p is to consider the largest minor, 
i.e., the determinant of p. Let me first examine how its 
sign depends on the eigenvalues of p. If conditions (J3]), (jl]) 
and ([B]) are satisfied, then by the interlacing inequality, 
p cannot have more than one negative eigenvalue. Thus, 
if det(p) is positive, then p is positive too. On the other 
hand, the negativity of det(p) implies that exactly one 
eigenvalue of p is negative. Finally, if det(p) = the 
interlacing inequality and condition ([3]) states that either 



Ai < 0, A 2 = 0, A 3 > 1, 



Ai = 0, < A 2 < A 3 < 1. 



(7) 



(8) 



However, in the first case the trace of the square of p 
should be greater or equal to the trace of p 

Tr(p 2 ) > Tr(p). 

One can easily check that 

3 



Tr(p 2 ) = £ 



5 £(4 +4)) < 



£ ( ^ ? + 2 £(<^ 

3>i 



= (£> =1, (9) 



therefore only the second case ((U can hold. Thus p is 
positive semidefinite if and only if all of its principal mi- 
nors are nonnegative. This ends the proof. 

To examine the structure of det(p), and additional con- 
straints on u>i, and let me consider the three 
two-dimensional vectors = (aij,/3ij), which can be 
simply interpreted as projections of Bloch vectors onto 
XY plane. For this purpose it is convenient to introduce 
the relative length of v~tj 



and the angles between and axes 

,'/V 

ipij = arctan 1 



(10) 



(11) 



The relative length is valid only for nonzero u>i and Wj and 
is bounded by < d^ < 1. Using these new parameters, 
the determinant of p is given by 



det(p) = uji<jJ2<jJ3 1 



£4 

i,j>i 



2^12^13^23 COS$ 



where 



(12) 



(13) 



Since u>i are always nonnegative, one can only consider 
the inequality 



1 



£4 

i,j>i 



2^12^13^23 COS <E> > 0. 



(14) 



The relative lengths dij are related to the pureness of 
the state. For example, similarly to the qubit case, the 
state is completely mixed if all dij = and all are 
equal. On the other hand, if only one parameter dij = 
1, the corresponding 2x2 principal minor equals zero, 
implying det(p) = 0. The other two principal minors are 
positive, therefore p has to have two positive eigenvalues. 
In this case the qutrit is confined in a two-dimensional 
subspace, thus cannot be maximally mixed. Moreover, 
the left hand side of (Tj"4")) is less or equal zero, thus the 
inequality is satisfied only if $ is a multiple of 2n and 
the other two relative lengths are equal. Finally, if two 
coefficients dij = 1, the third one has to equal 1 too, 
and again, $ has to be a multiple of 2ir. In this case all 
2x2 principal minors are 0. Moreover, Tr(p 2 ) = Tr(p), 
which implies that the qutrit is in a pure state. Note, 
that the pure state of the qutrit is defined by only four 
independent parameters, which in this case are two ip^ 
and two u>i. 

The pureness of the qutrit is also dependent on $. 
The parameter $ can differ from k2ir only if all dij are 
less than one. Manipulation with relative phases ip^ can 
change the rank of p from three to two. From the point 
of view of vectors vlj the dependence on $ tells about 
allowed angles between them. The greater is the rela- 
tive length, the more dependent on $ are the vectors. 
Some values of $ are not allowed for specific dij due to 
(|14p . This is an important constraint, because it states 
that for some relative lengths the vectors cannot point 
in arbitrary directions. If at least one vector has relative 
length equal to one, then 



<^12 - <^13 + ^23 = 0. 



(15) 



The set of states for which all values of $ are valid is 
given by the inequality 



i,j>i 



1 - £ 4 - 2d 12 rf 13^23 > 0, 



(16) 



which is obtained by taking cos$ = — 1. It is also clear 
that the dependence on $ vanish when at least one rela- 
tive length is zero. 
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FIG. 1: The Bloch vector for a qubit subspace of the qutrit. 



Top 

— the two-dimensional representation with normalized vec- 
tors, whose length is given in terms of dij , and with explicitly 
written maximal lengths 2^/uJiuJJ. Bottom — a schematic 
picture of the three Bloch vectors ulj 



IV. THREE BLOCH VECTORS 

The full set of constraints on , and uii is given 
by Eqs. ©, © and fTi|). One is able to represent 
the state of qutrit via three real two-dimensional vectors 
ulj = {puj , fiij ) , only if one also gives the values of the 
maximal allowed ranges dij, or simply the parameters u>i. 
However, a three-dimensional Bloch vector description 
can be used. Let me introduce the three new operators 

These operators are not linearly independent, since 

C12 - C 13 + C 23 = 0. (18) 
The expectation value of Cij is equal to 

(dj) =b} i -Uj=')ij. (19) 

Due to the commutation relation (fT7|) the three three- 
dimensional vectors ulj — (otij , (3ij , jij ) are the Bloch 
vectors of qubit subspaces, with vectors v~lj being their 
projections onto AijBij planes. 
The length of ulj cannot exceed 

+/%+ 7 £- < 1 - w fc = Rij , (20) 

where is the radius of the ij-th Bloch sphere. This 
is an additional important constraint. Also note, 

lij — Rik — Rjk- (21) 

One of such three vectors is depicted in Fig. [TJ 

The relations (|L8ti20[) lead to yet another two interest- 
ing constraints on vectors ulj 

712 - 7i3 + 723 = 0, (22) 



and 

\ui2\ + \ui 3 \ + \ui 3 \ <2. (23) 

The three-dimensional vectors ulj are even more depen- 
dent than the two-dimensional ones (namely vlj ) , but this 
is not surprising, since the number of parameters describ- 
ing the state is constant. In fig. [2] two representations 
of the state (5, e , ^) T are given. The vectors of the 
two-dimensional representation are normalized, i.e., their 
length is given in terms of dij . 

It is convenient to group all of the constraints on the 
Bloch vectors and to list them in a compact form: 



Si,j>j Rij — 2, 


< Rij < I, 


(i) 


I U*ij | — Rij 7 




(ii) 


lij = Rik — Rjk 




(iii) 


Rij =|[2 + (-l 


) i+ Hk + i-iyijk] 


(iv) 


712 — 713 + 723 - 


= 0, 


(v) 


i-E^4 + 


2^12^13^23 COS$ > 0. 


(vi) 



All the previous constraints on aij, fiij and uii are con- 
tained in the six listed above. This set is somehow over- 
complete, since some constraints can be derived from the 
other ones, however I believe that for the clarity it is 
convenient not to exclude them from the list. 

It is obvious that the set of all Bloch vectors for the 
qutrit is much richer than the one for the qubit. That 
is why in the discussion below I will narrow down to 
examine only the basic properties of the representation. 
The detailed study of mixed states is not presented here. 

A. Pure states and mixed states 

A nice and important property of both 3D and 2D rep- 
resentations is that for pure states all three vectors lie on 
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the corresponding Bloch spheres (circles) . This resembles 
the well known property of the qubit. Moreover, since a 
pure state is uniquely defined by only four parameters, 
it is sufficient to describe it by only two vectors, say u\2 
and U13 (alternatively v\2 and V\ 3 together with u)\ and 
UJ2)- The third vector is fully dependent on the first and 
the second one due to (i-vi). 

On the other hand, the completely mixed qubit state 
is represented by a null vector, i.e. a center of the Bloch 
sphere. Similarly, in the case of the qutrit one can eas- 
ily check that in the 3D representation all three vectors 
can vanish only if all three spheres are of equal radius 
2/3. However, this corresponds exactly to the completely 
mixed state. As already mentioned, in the 2D represen- 
tation a degree of mixture is related to the relative length 
dij . By analogy, in the 3D picture it is related to a ratio 
\ulj\/Rij. However, this ratio is not a proper measure of 
mixedness since it does not include the dependence on 
the angles tfij. 

Another interesting resemblance to the qubit Bloch 
sphere picture is the fact that if all three vectors are lying 
along Cij axes, the corresponding states are classical, in 
the sense that their density matrices are diagonal. Such 
states correspond to classical mixtures if the computa- 
tional basis consist of classically allowed states. 

B. Orthogonality 

Two states of the qubit arc orthogonal if their corre- 
sponding Bloch vectors are anti-parallel. For a given pure 
state of the qubit there is only one uniquely determined 
orthogonal state. It is not the case for the qutrit, since 
for an arbitrary pure state there is an orthogonal two- 
dimensional subspace containing infinitely many states, 
which arc not necessarily pure. Of course, two qutrit 
states can be orthogonal if at least one of them is pure. 
One can find an orthogonality criteria by a direct evalua- 
tion of Ti{pp'} = 0, assuming one density matrix is rank 
1, which leads to 

= 2 + £ ■ v7/, (24) 

i i,j>i 

or equivalently 

= 2 - E R v R ii + E u ^ ■ <• ( 25 ) 

ij>i i>j>i 

In case both states are pure, the later relation becomes 
= 2 - ^ Wlj\\ul/\ (1 - cos%) , (26) 

i,j>i 

where % is the angle ulj makes with ulj '. Note that for 
two pure states confined in a two dimensional subspace 
the above criteria reduce to the qubit orthogonality cri- 
terion. For example, if uo 3 = uo' z = then R12 = R[ 2 = 1 
and the pairs of vectors U13 and u\ 3 together with U23 
and U23 are pointing in the same direction along Cij 
axes, thus are parallel. As a result cos #12 = — 1. 



C. Algorithm for valid Bloch vectors 

To summarize this section I will write down a direct 
algorithm for a construction of any three vectors repre- 
senting a valid quantum state: 

1. Choose a radius < Rij < 1 of an arbitrary 
Bloch sphere and choose an arbitrary vector ulj — 
(aij,Pij,jij) lying inside this sphere. (As a con- 
sequence, all remaining R^ and 7^ are fixed - 
constraints (i-iv)) 

2. Choose the lengths of the remaining two vectors, so 
that they are confined in the corresponding Bloch 
spheres. 

3. Choose angles between the projections vlj and the 
corresponding Aij axes, so that (vi) is satisfied. 

V. DYNAMICS 

In the following section I am going to consider the uni- 
tary evolution of the qutrit from the point of view of the 
three Bloch vectors ulj. The dynamics of the qutrit is 
given by eight SU(3) generators, which in our case can 
be taken to be A^, Bij and two of the three Cij opera- 
tors. Any unitary operation on the qutrit can be decom- 
posed into a product of unitary operations Ug{9) — e lGe , 
where G is one of eight generators and 9 is a real param- 
eter. Therefore, it is enough to examine actions of these 
operators on the three Bloch vectors. To do this, it is 
convenient to introduce a table of commutation rules for 
the above operators (see Table 1). 

Let me start with the operators Cij, which generate 
rotations of the vectors ulj about Cij axes by the angle 
20. In addition, they also rotate vector ulk about Cik 
and ujk about Cjk by the angle 9. However, for C12 and 
C23, the rotation of U23 and U12, respectively, is done in 
the opposite direction (by —6). 

In addition to the three standard cyclic commuting 
triples {Aij,Bij,Cij}, for which the commutator of the 
two subsequent operators gives the third operator mul- 
tiplied by 2i, there are four extra cyclic commuting 
triples: {A 12 , A 13 , B 23 }, {A 12 , A 23 , B 13 }, {A 13 , A 23 , B 12 \ 
and {B12, Bi 3 , B 23 }, for which the commutator of the two 
subsequent operators gives the third operator multiplied 
only by i. Due to the cyclic commutation relation, the 
mean values of the operators in each triple can be consid- 
ered as coordinates of yet another Bloch-like vector, for 
which the operators A^ and B^ are generators of simple 
rotations. 

From the point of view of the already introduced Bloch 
vectors ulj, the generators A^ and B^ cause an oscilla- 
tion of the coordinates of ulk and ujk- However, a unitary 
operation cannot change a pure state into a mixed state, 
and vice versa, therefore, oscillations of coordinates of 
two vectors have to occur together with simultaneous os- 
cillations of the Bloch sphere radii Rik and Rjk- This 
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will make pure states having their Bloch vectors always 
on the surface of the corresponding Bloch spheres. In- 
deed, Aij generates the following change of the Bloch 
sphere radii 

Rik(fi) = \ {l+uJ k +A l3 sm(2e + 4>)), (27) 
RjkiO) = ~ {l + u k -A ij sva{2e + <j>)). (28) 
The action of £?,-., is similar 



Rik(0) = -(l+w fc + B«cos(2fl + 0)), (29) 
R jk (0) = ~ (1 + w fc - By cos(20 + <£)). (30) 



In the above, A; 



is a phase factor. 

The unitary evolution of the qutrit is given by a rota- 
tion of the Bloch vectors ulj and a resizing of the Bloch 
spheres. It is more complex than the simple evolution of 
the qubit, however in this case the number of generators 
is almost three times larger than in the qubit case. In 
the next section I will show that these generators have a 
natural interpretation for spin 1. 



VI. SPIN 1 



Spin 1 matrices S x , S y , and S z are cyclic commut- 
ing operators [Si,Sj] — iEijkSk whose eigenvalues are 
— 1, and 1. The first resemblance to the operators A^ 
and Bij is that both, spin 1 matrices and A^ and B«, 
have the same spectra. Moreover, the same cyclic com- 
mutation relations of some triples of A^ and B^j ma- 
trices were presented in the previous section. Indeed, 
if one represents spin 1 matrices in a basis spanned by 
their eigenvectors corresponding to the eigenvalue zero 
(the three vectors |5j = 0), i = x,y,z, are orthogonal), 
one obtains A 23l B\ 2 and A\ 3 respectively. In this case 
O i = \S i = 0){S i = 0\. 

Next, note that squares of the operators of the ij-th 
Bloch sphere obey the following relation 



At 



Bi 



cl 



I -O k 



(31) 



As a consequence, (Bf 2 ) + (A\ 3 ) + (A% 3 ) = R 12 + R13 + 
R 23 = 2. On the other hand, the well known angular 
momentum formula yields that (Sf.) + (Sy) + (Sf) = s(s + 
1), where s is the spin number, which for spin 1 gives the 
right hand side equal to 2. The above seems to lead to the 
following interpretation of the three Bloch spheres: each 
Bloch sphere describes a different spatial coordinate of 
spin 1 and the sphere radius corresponds to (Sf). The 
question is, what is the meaning of the other two Bloch 
sphere coordinates orthogonal to the one denoting S$? 



The operators CV, = O,; — Oj, therefore CV, — \Si = 
) = 0) (Sj =0\ = S]-S?. A difference of 



squares of two spin matrices is a generator of a certain 
type of spin squeezing known as two-axes countertwist- 
ing [llj]. In fact, the remaining three operators A\ 2 , B13 
and B 23 are also generators of two-axes countertwisting. 
Note, that they can be obtained from the first triple as 
anti-commutatators {£?i2,^4i3} = B12A13 + A13B12 = 
-B 23 , {#12,^23} = B 13 and {A 13 ,A 23 } = A 12 . To 
see that anti-commutators of spin matrices are two-axes 
countertwisting generators, let me write them explicitly 



{Si, Sj} — SiSj 
S, 



SjSi — 



$ i ~\~ S j Si ~\~ S j 

V2 V2 



Sj Si 



s 



V2 V2 



3 q2 

°i+3 



s 



0)<5i = 0| 



where Si+j and Si-j are spin matrices for the orthogonal 
directions i + j and i — j . 

The physical meaning of coordinates of a single Bloch 
vector for spin 1 are mean values of the following observ- 
ables: spin coordinate Si, generator of two-axes counter- 
twisting along two orthogonal axes j and k (which are 
both orthogonal to i), and generator of two-axes coun- 
tertwisting along axes j + k and j — k. At this point, let 
me discuss what information one gains while measuring 
the last two observables. The eigenvectors of a two-axes 
countertwisting generator Sf — Sj are \Si = 0), \Sj = 0) 
and \Sk — 0), where directions i, j and k are mutu- 
ally orthogonal. If the measurement yields an outcome 
\Si = 0), one knows that the mean spin vector is defi- 
nitely not lying along i axis, therefore the measurement 
of a two-axes countertwisting generator corresponds to 
the question: along which one of the three mutually or- 
thogonal directions i, j and k the spin vector is definitely 
not lying (see (l2l|). 

Finally, one also sees that every unitary transformation 
of spin 1 can be represented as a combination of three 
rotations and six two-axes countertwistings (however, it 
is enough to consider only five of them, since Cjj are 
linearly dependent). Moreover, an intuitive interpreta- 
tion of rotations and squeezing makes the change of the 
Bloch sphere radii quite natural, since rotations inter- 
change spatial coordinates and squeezing deforms them. 

Similarly to the qubit case, the proposed Bloch vec- 
tor representation of the qutrit has natural interpretation 
within the framework of the spin. It is an open question 
whether one can obtain similar property for the Bloch 
vector representation of higher dimensional systems. 



VII. HIGHER DIMENSIONS 

An analogical operator basis to the one used in this 
work can be easily constructed for any d-dimensional sys- 
tem. In such a basis one can easily represent any state 
by a collection of Bloch vectors. One may also try to 
find an analogy between basis operators and spin oper- 
ators, however the basis operators are either rank one 
projectors or Pauli matrices acting in a two dimensional 



7 



subspace, thus have eigenvalues 1, —1 and 0, which is 
d — 2 times degenerated. On the other hand, spin opera- 
tors, whose eigenvalues are 



2d— 1 



2(2-3 



2<2-l 



are 



2 , 2 , . . . , 2 

not degenerated and act on the whole space. Yet, they 
can be written as a direct sum of the operators propor- 
tional to the Pauli matrices. In case d is odd, one also 
has to include in the direct sum expansion, just like it 
is done for the qutrit. 

It is much harder to show what constraints are to 
be satisfied by the collection of Bloch vectors to prop- 
erly represent quantum states. It is obvious that the 
interlacing inequality requires all principal sub-matrices 
of the density matrix to be positive semidefinite. This 
implies that all principal minors of p are nonnegative. 
Unfortunately, for general density matrices the converse 
may not hold. Nevertheless, the non-negativity of prin- 
cipal minors gives necessary conditions for the positive 
semi-definiteness of p. Therefore, there are at least 
Y^rf ( 



2 d constraints on the basis expansion co- 



efficients, where the case n = corresponds to the nor- 
malization constraint. The exponential number of con- 
straints is in stark contrast to the quadratic number of 
qubit subspaces. Due to this fact, the higher is the level of 
the system, the more dependent the corresponding Bloch 
vectors are. There are d 2 — 1 parameters characterizing 
the state of the system and the average number of pa- 
rameters per Bloch vector is ■ 

The multi- vector representation can be especially use- 



ful for pure states, which are described by 2(d — 1) pa- 
rameters. In such case, the average number of parame- 
ters per Bloch vector is j an d one may expect the state 
to be representable by fever vectors. In particular, if it 
was possible to encode two parameters per vector, d — 1 
vectors would be enough to uniquely represent the state. 



VIII. CONCLUSIONS 

I have introduced a Bloch-like representation of the 
qutrit, which uses three vectors instead of one. The full 
set of constraints on vectors has been derived (Eqs. i-vi) 
and it was shown that many properties of the represen- 
tation resemble the ones for the qubit. In particular, 
pure states correspond to vectors lying on the surface of 
Bloch spheres and mixed states to vectors lying inside. 
Moreover, the dynamics of the qutrit can be described by 
rotations of the Bloch vectors and a change of the Bloch 
sphere radii. The generators of these transformations 
have natural interpretation for spin 1, namely, they gen- 
erate rotations and two-axes countertwisting spin squeez- 
ing. I also discuss basic properties of a generalization of 
this representation to higher dimensions. 
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«Bl2 


iA 13 
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A 23 
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iB 23 
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2iC*23 


-2iB 23 
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M13 
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—iA 12 
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2iA 23 


C23 


-iB 12 


iA 12 





iB 13 
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TABLE I: Commutation rules of Aij, Bij and dj. 



